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Abstract. The transport of reactive species brought about by ions traversing tissue-like medium is analysed 
analytically. Secondary electrons ejected by ions are capable of ionizing other molecules; the transport of 
these generations of electrons is studied using the random walk approximation until these electrons remain 
ballistic. Then, the distribution of solvated electrons produced as a result of interaction of low-energy 
electrons with water molecules is obtained. The radial distribution of energy loss by ions and secondary 
electrons to the medium yields the initial radial dose distribution, which can be used as initial conditions 
for the predicted shock waves. The formation, diffusion, and chemical evolution of hydroxyl radicals in 
liquid water are studied as well. 
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1 Introduction 

It is widely accepted that secondary particles such as elec¬ 
trons, solvated electrons, and free radicals play the key 
role in the radiation damage with ions P2II5] . These par¬ 
ticles are produced as a result of ionization and excitation 
of medium by ions as well as inelastic interactions of sec¬ 
ondary electrons of the first and following generations. The 
entire scenario of radiation damage with ions, considered 
in the multiscale approach m, includes the chemical 
stage when the above agents, formed in the vicinity of the 
ion’s path, interact with nearby biomolecules. Then fol¬ 
lows the predicted thermomechanical stage, during which 
the shock wave, generated by the gradient of energy de¬ 
posited in the medium, may rupture the covalent bonds 
of biomolecules, and transport the reactive species formed 
near the ion’s path to a much wider region. Then these 
species may interact with more distant biomolecules [T]. 

The analysis of transport of secondary particles is im¬ 
portant for all stages of this scenario. While the molecular 
dynamic simulation is the most adequate method of study 
for the latter stages [5] , an analytical random walk method 
can be applied for the first stage pQ. In this paper, we elab¬ 
orate the application of this method to describe the trans¬ 
port of two generations of electrons and reactive species 
such as free radicals and solvated electrons. We also anal¬ 
yse the location of energy deposition by secondary elec¬ 
trons and reconstruct the radial dose distribution based 
on the diffusion scenario. 


2 Transport of secondary particles originating 
on the ion’s path 


A large number of secondary electrons and radicals is pro¬ 
duced in the immediate vicinity of the ion’s path. The sec¬ 
ondary electrons are ejected as a result of ionizations of 
the molecules of the medium by projectiles. The radicals 
originate as a result of either molecule excitation by ions 
or as a result of hydrolysis due to local heating as a result 
of ion’s passage. The production of secondary particles on 
the ion’s path is described by the number of produced par¬ 
ticles of hue “j” per unit length of the ion’s path, dN p j/d(. 
An additional subscript “p” indicates that these particles 
originate from the ion’s path. 

The tree-dimensional diffusion of secondary particles 
such as electrons and free radicals is described by the fol¬ 
lowing equation, 


dn pj (r,t) 

dt 


DjX7 2 n pj (r,t) 


n p j(r, t) 


(1) 


where n P j( r, t) is the average number density of corre¬ 
sponding secondary particles at a location r (this vector 
connects points of origin of particles on the path and their 
observation), Dj is the corresponding diffusion coefficient 
in the medium, tj is the average lifetime. For electrons, 
Dj = D e = vl/ 6, where v is the electron’s velocity and l 
is their elastic mean free path in the medium. 

Equation HD has a well-known solution m 

dn P j(t 7 r) dN P j ( 1 \ 3 ^ 2 ( r 2 t \ 
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where r is the length of r. Equation m can be integrated 
over £ to obtain the number density for each species orig¬ 
inating at the ion’s path, 


n pj(t, p) 


dNpj 1 / p 2 t\ 

dC, AnDjt ^ y ADjt tj ) 


(3) 


The flux of secondary particles, d<P p j, through an ele¬ 
ment of a surface dA is given by 

d$Pi = {dA • VD 3 dn P^ r ^ dC (4 ) 


where the integration is done only over £, but not dA. In 
order to calculate the flux of secondary particles through 
a cylindrical shell of radius p and length £, with the ion’s 
path as the axis, dA is chosen to be —pd<f>£, where 4> is the 
azimuthal angle. The gradient of spherically symmetric 
dn pid’ r l i s collinear with r. Then, substituting the dot 
product r ■ p = rpp/r = p 2 in (pT|> . we obtain 


molecules are likely to become sources of free radicals. The 
latter also originate away from the path. After these ion¬ 
izations, both the newly ejected electrons and those that 
ionize the molecules have lower energies and, therefore, 
shorter elastic mean free paths. Finally, radicals propagat¬ 
ing from the path can react with the medium and produce 
new reactive species. All these effects can be described us¬ 
ing the same approach, described in this section. 

Each of these processes can be described by the cou¬ 
pled transport equations, 


dn pj (r,t) 


dt 


— DjS7 n pj { r,, 


drii(r,t) 

dt 


*)-£ 

= AV 2 »,(r,f) + ^ 


n pj{ r, t) 


‘ji 


(r ,t) m(r,t ) 


( 8 ) 


where index "i" marks the a specific kind of newly formed 
reactive species. The first of equations © describes the 
random walk of secondary particles ejected by the ion. 
This equation can be solved the same way as m to give 


dd'pj — 


r or aC 


dC (5) 


n e {t, p) 


dN e 1 / p 2 _ t_\ 

d £ 4-7T D e t ^ \ 4 D e t t J 


(9) 


Since this expression is cylindrically symmetric, there is 
no dependence of </>, and A can be taken as the area of a 
cylindrical shell of radius p and length £. Then d> p j is the 
total flux of radicals through this shell. 



Since the characteristic spatial scale in radial direction is 
in nanometers and in the axial direction is micrometers, 
d ^ J can be assumed to be constant. Then © can be 
integrated to give the flux, 


where 1/r = C V 7 *- 

Then, the rest of equations © become diffusion equa¬ 
tions with sources, and they can be solved analytically as 
well with a use of the Green’s functions, which are given 

by, 


Gi(t — r — r'> — _ g,)) 


3/2 


x exp ( — 
The solutions are: 


47t Di(t — t') 

(r-r ') 2 t-t' 


4 Di(t-t') r, 

it, r) = J G(t — t',r — r') 




( 10 ) 


( 11 ) 


‘Ppj 


dN pj (p,t) 

dt 


dN P j p 2 £ ( p 2 

-1C AD^ exp {-W~t 



which can be further integrated over time to obtain the 
total number of secondary particles of hue j, produced 
on the path, incident on the cylindrical shell. Such a cal¬ 
culation (in this work, only shown for this example) is 
important for the finding of the probability of a lesion in 
biomolecules [J. 


3 Transport of particles produced as a result 
of interactions of secondary particles with the 
medium 

Many of secondary electrons ejected from the location of 
the ion’s path have the energy sufficient for ionization of 
molecules of the medium. These ionizations take place 
at different distances from the path. Then, the ionized 


Each of these solutions requires four integrations, which 
are rather bulky, but doable. In order to simplify these 
integrations, it is advisable to use cylindrical coordinates 
with C-axis along the ion’s path. The expression looks as, 


r) 


! dNr 


3/2 


Tij d( J \AnDi(t-t') 

(r-r ') 2 t-t' 


x exp I — 


1 


4:irD e t' 


exp 


ADi{t-t') 

p' 2 t' 


AD P t’ 


Ti 

dt 1 dr' 


( 12 ) 


where 


(r _ r ') 2 = p 2 - 2 pp' cos(</> - 1) + p' 2 + (C - CO 2 - (13) 

and dr' = p'dp'dxj)'dC,'. Since dependence only appears 
in the exponential of (fT21) after substituting (fT51) . this inte¬ 
gration in infinite limits can be done first; and it is equal 
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to \J 4tt Dj(t — t'). Similarly, the dependence on the az¬ 
imuthal angle also only appears in the exponential of m 
after substituting m and 


/ 


exp 


/ 2p^_cos(^--^)\ 

V 4A(t - t') ) 


dip' = 2 ttIq 


PP 


2 Di(t - f) 



where Iq is a Bessel function. After this, two more integrals 
are remaining, 


rii(t, r) = 


1 dN e 


Tij dC, J 8nDiD e (t — t')t' 


x exp I — 


P 2 + P' 2 


t - t' p' 2 t' 


4A(t -t') Ti AD e t' Ti - 


PP 


2Di(t — t') 
First, the integral over p' is of the kind 


dt' p'dp'. (15) 


exp - 


P 


,12 


p'dp' = — exp 


(A ') 

\4& 2 J 


(16) 


After substituting the expressions for a and b into (flTH) . 
(usd becomes 


x exp 


rii{t,p ) 


1 dN e /■* 1 

47rr ij d( J 0 D e t' + Di(t — t') 


A{D e if + Di{t-V)) 


t — t' t' , . 

- dt'. (17) 


The last integration is not analytic, but it can be readily 
done numerically. Thus, nt is obtained as a function of 
time and the distance from the axis; it can be added to 
n P i for a given species, given by ([3|). 

The contribution to the flux due to Vn^ can also be 
calculated as a function of p and t. All these calculations 
include parameters dN p i/d(, Ti, and for each species 
of interest. 


4 Transport of electrons ejected by secondary 
electrons 


Many of the secondary electrons ejected have enough en¬ 
ergy to ionize water molecules in the medium. Indeed, 
more than 65% of these electrons have energies higher 
than the ionization potential of water molecules. The sec¬ 
ond “wave” of ionization has been discussed [7| , but only 
from the point of view of remaining energy, and it was 
concluded that the third wave is insignificant. For some 
applications, it is sufficient to consider the diffusion of sec¬ 
ondary electrons leaving the second wave aside [8]. How¬ 
ever, Sec. [3] suggests an analytic method for accounting of 
the second generation of electrons. 


The formation and transport of the second wave ion¬ 
ization can be described by equations ©, 


dt 


= £>iV 2 m(r,f) - 




Tl 


dn 2 (r,t) 


dt 


= -D 2 V 2 n 2 (r, t) 


ni(r,t) n 2 (r,t) 


Tl 


T 2 


(18) 


Here, index “1” marks secondary electrons of the first gen¬ 
eration. Their energy is taken to be equal to be the average 
energy of secondary electrons formed in the vicinity of the 
Bragg peak, i.e., 45 eV. Hence, D i =vl /6 = 0.265nm 2 fs _1 
and Ti = hon/v = 0.64fs (all mean free path data are 
taken from Refs. mm- Index “2” corresponds to 15-eV 
electrons, formed as a result of the second wave of ion¬ 
ization by secondary electrons; D 2 = 0.057nm 2 fs _1 and 
r 2 = 15.3fs. The solutions of the above equations are ob¬ 
tained similarly to those of ©, 


ni(t,p) = 


dNe 


dC, AnDit 


exp - 


n 2 (t,p) = 2 


ADit 
1 


t 

Tl 


x exp 


1 dN e r 

47TT1 d( J 0 Dit' + D 2 (t-t') 
p 2 t-U 


A{Dit' + D 2 (t-t')) r 2 


- ^ dt'. (19) 


Notice, that n 2 (t, p) is doubled since it comprises the newly 
ejected electrons as well as the secondary electrons of the 
first generation that lost energy. The results for these num¬ 
ber densities are shown in Fig. |T| From this figure, it is 
clear that after a few fs after the ion’s passage, all sec¬ 
ondary electrons lose energy and the number density of 
secondary electrons is by and large is given by n 2 (p,t). In 
Fig. El n 2 ( p , t) is plotted at different times as a function 
of distance from the path. With time the distribution be¬ 
comes a little broader, but the main effect is the exponen¬ 
tial decrease with time. As a result of these decrease, the 
so-called pre-solvated electrons are formed. A pre-solvated 
stage of electrons is a transition stage between low-energy 
ballistic electrons and a relatively stable compound of elec¬ 
trons with water molecules known as solvated electrons. 
This transition takes about 1 ps. 

The formation of pre-solvated electrons can be found 
from the conservation of electrons corresponding to m 


dn aq (r,t) 

dt 


n 2 (r,t) 

t 2 


( 20 ) 


where the subscript “a<f corresponds to pre-solvated elec¬ 
trons. In (1201) . the diffusion and all chemical terms are 
dropped since the diffusion ( D aq = 4.5x 10 _6 nm 2 fs _1 jTTj ) 
and chemical reactions are happening on a much longer 
scale. Equation (I2U1) can be integrated over time to get the 
initial distribution of pre-solvated electrons. During the 
time of formation of solvated electrons, pre-solvated elec¬ 
trons may interact with nearby biomolecules, but, since 
their diffusion coefficient is much smaller compared to 
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t(fs) 

Fig. 1. (Colour online) Number densities of secondary elec¬ 
trons of second generation at 0.5- (solid) 1-, and 1.5-nm dis¬ 
tances (dashed with decreasing dash size, respectively) from 
the path as functions of time. The distributions of the first 
generation at the same distances are shown in the inset. 



p(nm) 


Fig. 2. (Colour online) Number densities of secondary elec¬ 
trons after they have ionized water molecules at 5, 10, 15, and 
20-fs times (solid and dashed lines with a diminishing dash 
size, correspondingly) as functions of the distance from the 
ion’s path. 



P(nm) 


Fig. 3. (Colour online) Number density of pre-solvated elec¬ 
trons (solid line) and initial distribution of hydroxyl radicals 
(dashed line) as a function of distance from the ion’s path at 
50 fs when the transport of ballistic electrons is over. The ratio 
of chemical annihilation of hydroxyl radicals to their diffusion 
as a function of the distance from the ion’s path is shown in 
the inset. 


those of ballistic electrons, their location does not change 
considerably compared to the position of last inelastic 
event. Therefore, the distribution obtained by the inte¬ 
gration of (BUI) is the initial radial distribution of solvated 
electrons. This distribution is shown in Fig. [3] Later on 
they slowly diffuse and react with other solvated electrons 
and hydroxyl radicals to form stable OIT - nu. Chemical 
reactions dominate the diffusion and, as will be discussed 
below, given that diffusion is the only mechanism for the 
transport of solvated electrons, they are unlikely to leave 
the track. 

It is also possible to assess the energy deposition den¬ 
sity (dose), e(p,t), in the medium by ions and secondary 
electrons. In order to do this, we assume that the average 
energy, w, is deposited to the medium with each ioniza¬ 
tion. Then, the rate of energy deposition is proportional 
to the rate of inelastic events, 

= m ( wJ( () + HiAi! + 5M\ 21) 
at \ aQ ti 72 / 

The first energy deposition occurs right at the ion’s path 
where the molecules are ionized by the ion; this corre¬ 
sponds to the first term on the right hand side of m 
where (5’s are the corresponding (5-functions. The second 
deposition (second term on the r.h.s.) is the ionization by 
secondary electrons at the end of their ionization mean 
free paths. Finally, the third deposition is due to remain¬ 
ing energy loss due to excitation of molecules by the dec- 
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Fig. 4. (Colour online) The radial dose (in MGy) distributions 
as functions of the distance from the ion’s path at times 50, 20, 
10, and 5 fs are shown with solid and dashed lines with a dimin¬ 
ishing dash size, correspondingly. The solid line also represents 
the initial distribution of hydroxyl radicals. The corresponding 
labels are shown on the right side of the frame. The dots mark 
the radial dose data for 2-MeV/u }l2j multiplied by the factor 
of four. 

trons of second generation. After that, the electrons enter 
the pre-solvated stage. 

The time integration of (EH) gives the dependence of 
the radial dose on time. The radial dose distributions at 
times 5, 10, 20, and 50 fs are shown in Fig. [I] At small 
radii, the distribution is due to primary ionization and it 
does not change with time. At larger radii, the dose slowly 
increases because of energy loss by the second generation 
of electrons. The shown results are obtained with w = 
16.5 eV, which corresponds to the normalization 

J j ^ dt2n P d P = Se ’ ( 22 ) 

where S e = 900eV/nm is the linear energy transfer (LET) 
of a single ion at the Bragg peak for carbon ions [7]. In 
Fig. 01 the results of integration of (ITU are compared with 
the radial dose of Ref. [12]. While the shapes of these dis¬ 
tributions are alike, the absolute values are different by 
the factor of about four. Part of this disagreement (fac¬ 
tor of 1.9) can be explained by the fact that in [T2], the 
radial dose presented for 2-MeV/u, i.e., proximal of the 
Bragg peak for carbon ions. The remaining factor can be 
due to the energy straggling present in Ref. m , but ab¬ 
sent for a single ion data of EH- Also, while the solid line 
represents the radial dose at 50 fs, in Ref. [T2] there is 
no information about the time. This is typical for Monte 
Carlo simulations, but since they are compared with ex¬ 


perimental data the time corresponding to dots is likely 
to be on a ps scale. 

The results shown in Fig. 0] may be used as a starting 
point for shock wave development mm, since the radial 
dose distribution evolved by the time < 50 fs determines 
the temperature distribution around the ion’s path fol¬ 
lowing from the process of energy relaxation, in which 
the energy stored in electronic excitations is transferred 
into vibrational excitations and then to translational de¬ 
grees of freedom. The characteristic time for the decay of 
an electronic excitation in Na clusters is estimated to be 
about 0.4 ps [T3] and it is likely to be several times longer 
for the liquid water, which is consistent with the analysis 
performed in Ref. [141 . The dependence shown in Fig. [4] 
is somewhat similar to that discussed in Ref. m, but the 
difference is that EH includes diffusion of secondary elec¬ 
trons, but does not include thermal conductivity, while 
Ref. P3] does the opposite. 

Finally, the radial dose has been studied in Ref. [83 
using a random walk approximation. In that work, we 
obtained a reasonably good comparison of shapes with 
Ref. [12 for intermediate distances. Only one generation of 
secondary electrons was used, but the parameters, such as 
mean free path and the relaxation time were chosen close 
to those for the second generation of this work. Therefore, 
the results obtained in Ref. [8] are meaningful, however, 
they do not include term that comes from 

primary ionization events. 

5 Transport and production of hydroxyl 
radicals 

Hydroxyl among other radicals plays the most important 
role in the DNA damage mm- Therefore, we are consid¬ 
ering the production and transport of these radicals in 
more detail. In the liquid water environment irradiated 
with ions such as carbon, the hydroxyl radicals are pri¬ 
marily produced as a consequence of ionization of water 
molecules, 

H 2 0 + + H 2 0^H 3 0+ + 0H', (23) 

and by excited water molecules, 

H 2 0*^H+OH. (24) 

Both of these processes can be initiated by either ions or 
secondary electrons that ionize and excite water molecules. 
The time scales play an important role in further analy¬ 
sis. The ionization is happening in 10 -17 s and secondary 
electrons diffuse by several nm within several fs. How¬ 
ever, OH' is produced through (E3l) and (HHl) only on a 
ps scale, i.e., after the electron transport is over and pre- 
solvated electrons are formed. Neither H 2 0 + , nor H 2 0* 
are remaining at rest during this time; they are diffusing 
away from the path. Since their diffusion coefficients are 
of the order of 10~ 6 nm 2 fs _1 , they diffuse by only about 
a/ DT form « Vio- 6 io 3 = 3 x nr 2 nm away from places 
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where the inelastic event (ionization or excitation of wa¬ 
ter molecule) has taken place. This distance is too small, 
therefore, we can consider the loci of inelastic events to be 
the initial loci of hydroxyl radicals. 

The distribution of in elastic events was analysed in 
Sec. 0] in Eq. CD, where the rate of energy deposition 
was obtained. In the first approximation, we can assume 
that each inelastic event produces a hydroxyl radical, then 
the equation, similar to ED for production of radicals is 
as follows, 

dn OR (p,t) = dN e j( 2 ) (pWf) + ni (P'V + n 2(/M) , 25 j 
dt d( T\ r 2 

Similarly to ED^ it can be integrated over time to obtain 
the initial number density of hydroxyl radicals, uoh(p)- It 
is shown with a solid line in Fig. 0] with labels on the right 
side of the graph. The results of the integration of second 
and third terms are shown with a dashed line in Fig. 0] 
to compare the distributions of hydroxyl with solvated 
electrons at larger distances. 

Thus, with the above assumptions, each ionization on 
the ion’s path produces one 45-eV electron, which after an¬ 
other ionization loses energy and produces two 15-eV elec¬ 
trons, which then become two solvated electrons. These 
two ionization events together with final energy loss by 
two 15-eV electrons produce four hydroxyl radicals. How¬ 
ever, the initial radial distributions of solvated electrons 
and hydroxyl are different, the hydroxyl strongly domi¬ 
nates near the path and they are the same at p > 1.5nm. 
Now, let us analyse the further evolution of hydroxyl. 

As can be seen from Fig. 0J the ballistic electrons are 
gone after about 50 fs, a time much shorter than that of 
formation of hydroxyl. Thus, by the time hydroxyl radicals 
are formed, there are no more active sources since both 
ionization and excitation processes are over. All reactions 
that eradicate hydroxyl are second order reactions with 
respect to perturbations in water; these are 

OH' +H H 2 0, 

OH + e aq OH , 

OH' + OH' H 2 0 2 , (26) 

and others m- This means that all equations such as m 
and E| acquire additional terms on the right hand side 
that account for chemical reactions, but lose terms like 
nou(j,t)/ tou, since hydroxyl does not react with water 
molecules and is stable by itself. For example, Eq. (0| be¬ 
comes 

— ^ = -DoHV 2 noH(r, t) — ^ fcomnoH^i, (27) 

where kom is the rate constant corresponding to the re¬ 
action of species with hydroxyl. This is a nonlinear 
equation, and we start with comparing all terms on the 
right hand side of (1271) . 

First, let us consider the hydroxyl radicals that are 
formed due to ionizations or excitations by ions. Since 
the diffusion coefficient for hydroxyl is so small (2.8 x 


10~ 6 nm 2 fs -1 ) compared to electrons discussed above, the 
first term in Eq. (1271) (for carbon ions in the vicinity of 
Bragg peak) is of the order of 10 if we consider the esti¬ 
mated diffusion radius to be 0.03-nm. The second term 
is equal to &oh,oh^oh’ s i nce hydroxyl by itself is the 
dominant reagent in this region. With /coh.oh = 1.0 x 
10 _5 nm 3 fs _1 [TT1I3112]. this term exceeds 200nm _3 fs _1 , since 
the concentration near the path is so high. This concentra¬ 
tion decreases by the factor of 100 within 2 ps as most of 
the hydroxyl becomes involved in the formation of perox¬ 
ide. Therefore, the hydroxyl, formed by direct ion’s action 
by and large does not exit the a sub-nm region around the 
ion’s path. 

The situation with hydroxyl formed by secondary elec¬ 
trons is somewhat different, since it is formed with a much 
smaller concentration and the diffusion term may be com¬ 
parable with that of peroxide formation. The equation of 
interest can be written as, 

dn 0 n(p, t) 1 d ( dn 0 u{p,t)\ 

— m — =Dou - p y p dp ) 

—kou,e aq noH(p,t)n aq (p,t) — /coh,oh^oh(/°: *)■ (28) 

This equation includes two relevant chemical reactions 
from (1261) . In the beginning, at t « lps, the number den¬ 
sity of hydroxyl, noH(p, t), exceeds that of solvated elec¬ 
trons, n aq (p,t), as can be seen in Fig. 0] The reaction 
constant /coh e aq is larger than /coh oh by the factor of 
4.2 mm®- Therefore, the reaction terms in (1751) domi¬ 
nate the diffusion term in the whole domain of p. The most 
conservative estimate of the ratio of reaction to diffusion 
terms in (1751) is shown in the inset of Fig. 0] It is apparent 
that terms become comparable only at p > 2.5 nm. Only 
a small (less than 10%) fraction of hydroxyl diffuses to 
these distances. 

Therefore, we can summarise that in the first approxi¬ 
mation, which can be improved by introduction of proba¬ 
bilities of formation of hydroxyl radicals, the initial distri¬ 
bution of OH' is shown in Figs. 0J0J These radicals would 
react with biomolecules such as DNA if the latter appear 
to be in the track. However, if the collective transport 
with a shock wave does not take place, the radicals react 
with each other and with solvated electrons and by and 
large do not leave the track. Other species, such as H 2 0 2 
and OH - may propagate on larger distances and may be 
detected outside tracks. 


6 Conclusion 

The random walk approximation, applied to two gener¬ 
ations of secondary electrons showed that most of them 
do not spread beyond the 2-nm cylinder around the ion’s 
path. Pre-solvated electrons, important agents of strand 
breaks in DNA molecules, are also formed within this 
small region. Their small diffusion coefficients do not allow 
them to be transferred far enough since they readily react 
with hydroxyl radicals abundant in the same region. The 
latter also react with themselves and their high number 
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density makes this (nonlinear) reaction much more impor¬ 
tant than their (linear) diffusion. Thus, we can conclude 
that the diffusion mechanism does not allow most of the 
reactive species to leave the few-nm track radius. In this 
work, carbon ions near the Bragg peak were used as an 
example. The effect of containment of reactive species in 
the track will be stronger for heavier ions and weaker for 
lighter ones, since the concentration of reactive species is 
much smaller. 

Alternatively, another transport mechanism, i.e., col¬ 
lective transport in a shock wave should be studied both 
theoretically and experimentally in order to understand 
whether the reactive species actually propagate to larger 
distances. This question is very important for the assess¬ 
ment of radiation damage with ions [Tj. The radial dose 
distribution obtained in Fig. [I] gives the initial conditions 
for the development of cylindrical shock waves. 
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